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Abstract 

(N 

C^i . The Regge- Wheeler equation describes the axial perturbations of Schwarz- 

schild metric in linear approximation. We present its exact solutions in terms 
of the confluent Heun's functions, the basic properties of the general solu- 
tion, novel analytical approach and numerical techniques for study of different 
boundary problems which correspond to quasi-normal modes of black holes 
and other simple models of compact objects. We depict in more detail the 
exact solutions of Regge- Wheeler equation in the Schwarzschild black hole 
[ interior and on Kruscal-Szekeres manifold. 

(N 
O 

q ; 1 Introduction 

Q" 1 . The well known Regge- Wheeler equation (RWE) 

99- d 2 t $ s>l + + v s> i) $ M = o (i.i) 



>< 



describes the axial perturbations of Schwarzschild metric in linear approximation 



and plays an important role in modern perturbation treatment of Schwarzschild 



black hole (SBH) physics. Its study has a long history and significant achievements 
®- 

The effective potential in RWE (11 .11) reads 



K,(r) 



x _ 1\ ( 1(1 + 1) W 



The area radius r > can be expressed explicitly as a function of the "tortoise" 
coordinate x = r+r Sch In (\r/r Sck — 1|) using Lambert W-function: r = W(±e x ~ 1 ) + 1. 
In the last formula sign " + " stands for SBH exterior r G (1, oo) and sign " — " — for 
the SBH interior r e (0, 1). Hereafter we are using units in which the Schwarzschild 
radius r Sch = 2M = 1. The most important case s = 2 describes the gravitational 
waves and is our main subject here. 
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The ansatz <E> S) j(£, r) = R LJ , S) i(r)e' luJt brings us to the stationary problem in the 
outer domain r > 1: 

d?R w>s>l + (u 2 - V a ,i) K,s,i = 0. (1.2) 

Its exact solutions were described recently [3] in terms of the confluent Heun's func- 
tions HeunC(a, /3, 7, 5, 77, r) [I]. 

In the exterior domain the variable r plays the role of 3D-space coordinate and 
the variable t - the role of exterior time of a distant observer. Because of the 
change of the signs of the components gu = —^/g r r = 1 — 1/r of the metric, in 
the inner domain the former Schwarzschild-time variable t plays the role of radial 
space variable r in = t G (0, 00) and the area radius r G (1, 0), plays the role of time 
variable. As a result, the Regge- Wheeler "tortoise" coordinate x G (— 00, 0) presents 
a specific time coordinate in the inner domain. For the study of the solutions of 
( 11.11) in this domain it is useful to stretch the last interval to the standard one by a 
further change of the time- variable: x — * t in = x — 1/x G (—00, 00). 

These notes are important for the physical interpretation of the mathematical 
results. In particular, the natural form of the interior solutions of (11.11) is 

where 




The dependence of this solution on the interior radial variable r in is simple. Its 
dependence on the interior time t in is governed by the (11.21) with interior-time- 
dependent potential V s j. In despite of this unusual feature of the solutions in the 
SBH interior, this way we obtain a basis of functions, which are suitable for the 
study of the corresponding linear perturbations. 

2 Exact Solutions of Regge- Wheeler Equation in 
Terms of Confluent Heun's functions 



In r variable the RWE reads: 



dr 2 r(r — 1) dr 
u 2 r 4 - 1(1 + l)r 2 + (1(1 + 1) + s 2 - 1) r + 1 - s 2 



... ., ■ ., Rw,s,i 0. 

r z (r — 1 r 
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The anzatz -Roj,s,z( r ) = r s+1 (r — l) llJJ e lujr H (r) reduces it to the confluent Heun equa- 
tion [1]: 

d 2 H ( 8 + 1 7 + 1\ dH 
+ la + - + - — + 



r(r 



dr 2 \ r r — 1 J dr 

( (* + \a{(3 + 1 + 2))r + 77 + ^ + a) 09 + 1)) IT = (2.2) 



with the specific parameters: a = 2za;, /3 = 2s, 7 = 2za;, 5 = 2u; 2 , 77 = s 2 — /(/ + 1) and 
three singular points in the whole complex plane C r [3]. Two of them: r = and 
r = 1 are regular and can be treated on equal foods. The third one: r = 00 is an 
irregular singular point. Note that after all the horizon r = 1 turns to be a singular 
point in the above sense, nevertheless with respect to the algebraic invariants of 
the Riemann's curvature tensor it does not define a singular surface in the very 
Schwarzschild space-time manifold M^ 1,3 ^{py(t, r, 9, (f>)}. 

Correspondingly, one can write down three pairs of exact local solutions - one 
regular and one singular - around each singular point X = 0, 1, 00: 

1) Frobenius type of solutions around the points r = and r = 1: 

*%(t,r) = ^R%{r), (2.3a) 

f tt(0)+ ( \ 

R {0)± ( r ) = r s+l e iu>r+^Hl-r) J u>,s,lK h 

"' sA [H eunC \2iu, 2s, 2iu, 2u 2 , s 2 - 1(1 + 1), r) ; 

*%(t,r) = e^R^r), (2.4a) 

^jWi ^ __ r s+lgiwr±iwln(r-l) x 

HeunC {-2iu, ±2iu, 2s, -2uj 2 , 2uj 2 + s 2 - 1(1 + 1), 1 - r); (2.4b) 

2) Asymptotic Tome type of solutions around the point r = 00: 
^£(t,r)=e^R^(r), (2.5a) 
<5» ~ e^" 1 - £ of = 1. (2.5b) 



<' r 



These solutions describe ingoing and outgoing waves, associated with the corre- 
sponding singular point, according to the following scheme. 



The solutions own the following properties: 

lim|S<J+(t,r)| = 00, (2.6a) 
lmx\^-(t,r)\ = 0, (2.6b) 



Table 1: In-out properties of local solutions. 
Ingoing solutions: Outgoing solutions: 



(0)+ 

LO,S,l 

,(!)+ 



(t,r) : 
(t,r) : 
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with the important justification in the case |u;| > 0, as described by the equations 

(2.7a) 

(2.7b) 

and illustrated in Fig{TJ The general solution of the RWE equation (11.21) can be rep- 



lim l<i!(r) 



oo, if arg(r) + arg(cj) G (0,7r)/mod 2n; 

0, if arg(r) + arg(cj) G (— 7r, 0)/mod 27T. 

0, if arg(r) + arg(co>) G (0,7r)/mod 2n; 

oo, if arg(r) + arg(c<j) G (— 7r,0)/mod 27T, 
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Figure 1: Domains of different behavior of the functions -R^]f (r). 

resented in the form -R^ iSj /(r) = Cx+R *si ( r ) + Cx~R^i~i( r ) "with proper constants 
Cx±- The different local solutions are related according the formula 



<g(r) 



rf^(u,sJ)R^Vr)+T*±(u,s,l)R 



X±, 



The main obstacle for pure analytical treatment of different interesting problems, 
related to the Regge- Wheeler equation, is that the transition coefficients r y± (u;, s, I) 
for X, Y = 0, 1, oo; 1^7 are not known explicitly. 
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3 Different Boundary Problems for the Station- 
ary Regge- Wheeler Equation 



Due to the existence of three singular points of RWE, we have three different two- 
singular-ends boundary problems [I] in the complex plane C r , FigEJ These describe 
different physical problems, related with SBH exterior, SBH interior and Kruscal- 
Szekeres (KS) extension of the Schwarzschild solution [5]. 

Besides we can consider regular-singular-ends boundary problems on the interval 
(r*, oo), r* > 1 which can describe the perturbation of the metric around compact 
spherically symmetric matter objects in rest. 



r=l 




Figure 2: Different types of two-singular-ends boundary problems. 



3.1 Perturbations of the SBH Exterior 

Most well studied are the linear perturbations of the SBH exterior [H[2]. Since the 
quasi-normal modes (QNM) are defined as a stable solutions of RWE (11.11) which 
describe pure ingoing waves both into the horizon r = 1 and into the space infinity 
r = oo, we can describe them using the solutions R^i( r ) 03) of (OIL which 
obey the non-explicit equation r^ l(o;,s, I) = 0. It yields the QNM frequencies 
uj = oj(n, s, I), n = 0, 1, 2, . . . . Our exact approach gives the spectral condition [3]: 



lim 

|r|^+oo 



r s+1 HeunC (-2iu, 2tu } 2s, -2u 2 , 2uo 2 +s 2 -1(1 + 1), l-| r | e -^ +argM ) 



0. (3.1) 



It expresses in an explicit way the condition that the ingoing at the horizon waves 
do not contain coming from the infinity waves and can be solved numerically using 
the computer package Maple 10. 
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The equation ( 13. ip reproduces with a great precision the known numerical results, 
obtained by different methods [2]. For example, for the basic QNM-eigenvalue our 
calculations give u(0, 2, 2) = 0.747343368836 + i 0.177924631360 [3], thus justifying 
the previously published results [2]. The first six eigenvalues are shown in FigJHJ 

The amplitude of all QNM solutions R^^isi) s i( r ) increase infinitely around the 
horizon and near the space infinity, according to the (I2.6aj) . This shows that Regge- 
Wheeler perturbation theory is not applicable in the corresponding vicinities of these 
points. 

3.2 Perturbations of the SBH Interior 

The perturbations of SBH interior were considered qualitatively for the first time in 
[6] . We present here an explicit exact treatment of the problem. 

In the classical theory of BH it seems reasonable to exclude from the most of 
the considerations the inner domain, because it does not influence directly the ob- 
servable outer domain. Nevertheless, the excitations of degrees of freedom in the 
inner domain may be essential for the consideration of the total BH entropy and 
for construction of the quantum theory of BH. The interior excitations may yield 
observable effects in the exterior of the BH, due to the quantum effects [TJ. 

Another reason to study the BH interior was discovered in the recent article 
[8] : It was shown that without excision of the inner domain around the singularity 
one can improve dramatically the long-term stability of the numerical calculations. 
Although the region of the space-time that is causally disconnected can be ignored 
for signals and perturbations traveling at physical speeds, numerical signals, such 
as gauge waves or constraint violations, may travel at velocities larger then that of 
light and thus leave the physically disconnected region. 

The perturbative treatment is related with the metric's dynamics and corre- 
sponding initial value problem: using the field equation one has to determine the 
evolution of the initial data on a proper Cauchy surface. In the exterior domain so- 
lutions of (ll.ip . which grow up infinitely with respect to the future direction of the 
exterior time t do exist. They are unstable and physically unacceptable. Therefore 
such solutions are excluded from any physical consideration [2J [Tj . 

Then one has to formulate a proper boundary problem for (II. 2p . which has a 
reach variety of solutions and may describe different physical problems. For a correct 
formulation of a given physical problem one restricts the class of the solutions using 
proper boundary conditions. Since the physics at the boundary influences strongly 
the solutions of the problem, the boundary conditions are a necessary ingredient of 
the theory and define the physical problem under consideration. Only after that the 
physical problem is fixed and one can prove the stability of the remaining solutions 
(QNM) and small perturbations of general form with respect to the exterior-time 
evolution [9]. Our point is to consider different type of mathematical boundary 
problems, related with (11. 2ft . 

By analogy with the consideration in the outer domain of SBH, in the inner 
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one we are considering only solutions, which are stable with respect to the future 
direction of the interior time. Therefore we study the solutions of (ll.2p . which are 
regular at the point r = 0. These have a nontrivial spectrum with novel properties. 
Since we are working in the framework of perturbation theory, the presence of "bad" 
solutions Ruti( r )i see Q2.3bh . which diverge and even are not single- valued, shows 
that, in general, two different types of situation can take place: either the considered 
here perturbation theory does not work in some domain around the interior-future 
infinity, or the full nonlinear problem has "bad" solutions, which have to be ex- 
cluded from the physical considerations using proper additional conditions. It is 
impossible to resolve this problem in the framework of the linear perturbation the- 
ory. Nevertheless, it is clear, that the opposite case of the regular solutions R^ s j(r), 
see ( I2.3bl) . has always a well defined physical meaning around future infinity of the 
interior time and deserves a corresponding study. For this purpose we are using a 
second pair of solutions around the horizon r = 1: 

£«± (r)=r^ 1 e ±w(Hl ^ hr) HeunC(±2iu,±2iu;,2s,-2u 2 ,2u; 2 +s 2 -l(l + l),l-r), (3.2) 

which can be derived from the solutions (12.4bl) . applying proper transformations of 
the Heun's functions jl]. For u) ^ the functions R^i(r) (13. 2\i define a pair of 

linearly independent solutions with the symmetry property R^fi( r ) = R^-usi( r )- 

1. For real uj 2 > the problem has a real continuous spectrum. The corre- 
sponding solutions Rw,2,2(tin) are illustrated in Figj3] They define a basis of stable 
normal modes for perturbations of the metric in the SBH interior and can be used 
for Fourier expansion of perturbations of more general form. 
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Figure 3: Two solutions of the continuous spectrum. 

2. The zero eigenvalue uj = yields a degenerate case in which the two functions 
(13.21) coincide. Then we obtain the infinite series of polynomial solutions: 

R (O J: , s ,i( r ) = Pochl72s,l-s) rS+1JacMP{l ~ S > 2S ' °' 1_2r) = r3 + - • • ' (3 ' 3) 
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which are finite on the whole interval r G [0, 1]. These functions form an orthogonal 
basis of stable normal modes. Here JacobiP(n, a, b, z) stands for the standard Jacobi 
polynomials and Poch(c, d) - for the Pochhammer symbol. 

3. The discrete spectrum of the problem at hand corresponds to pure imaginary 
values of the parameter u = iujj 0,[3], i.e. to Laplace transform of the perturbations 
of general type. Thus, studying the the case u 2 < we are constructing a basis for 
Laplace expansion of the perturbations of more general form in the SBH interior. 

Then the functions Ru7i( r )> R^jfi( r ) anc ^ a ^ their parameters are real. In 
particular, 

< s » =r?>, s, l)R {1 JXi(r) + r?l(u;, s, l)R ( lUr) (3.4) 
~ cos a R (1 JXi(r) + sin a R^/r) 

with real transition coefficients r5^(co>, s, /) and real mixing angle a, defined via the 

relation tan a = r°:(cu, s, 0/r?;(cu, s ,i)eR for 

The physical meaning of the mixing angle a is obvious: since the solutions 
Rujti( r ) describe the ingoing into the horizon waves and the solutions R\J s i(r) - 
the outgoing from the horizon waves, this angle describes the ratio of the amplitudes 
of these waves in their mixture (13.41) . The value a = corresponds to absence of 
outgoing waves (i.e. to the authentic SBH boundary conditions at the horizon). The 
value a = 7i/2 describes the case without ingoing waves. The values a G (0,7r/2] 
describe an extension of the simple SBH problem in which the interior domain has 
a more general meaning then a pure SBH interior. Thus, fixing the value of the 
mixing angle a G [0, 7r/2] we are defining completely a more wide class of problems 
than the traditional SBH. 
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Figure 5: The first 18 
Figure 4: The inverse po- eigenvalues (for n=0, 

tential -V(r). 2,..., 17, s=l=2 and 

a = 0). 
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The discrete imaginary spectrum can be obtained from the spectral equation: 



in: 



,(o)- 

U),S,l \ 



cos a -R^7z( ri ) +sin a i? 



(i)- 

U!,S,l 
(1)- 



cos a R y J !S j{r 2 ) +sin a R>J ta>l {r 2 



(3.5) 



where r 12 € (0,1), 7^ r 2 are two otherwise arbitrary points. Because of the 
symmetry property of solutions R^ti( r )> ^ ^ s enou gh to study only the solutions 
u) n , s ,i of the equation (13.51) with positive uj = Im {uj n ,s,i) > 0. 

In the Figj5]we see the first 18 eigenvalues (including the zero eigenvalue) of the 
RWE for SBH interior for s = I = 2 and a = 0. The two series: n = 0, 1, ... 6; 
and n — 7, . . . , which correspond to the two wells of the inverse potential — V(r), 
are clearly seen. They are placed around the straight lines uj = n/2 + 1/2 and 
u)i — n/2 + 1, correspondingly. 

In Fig J6] we see a specific attraction and repulsion of the eigenvalues ui nySy i{a). 
Such behavior is typical for the eigenvalue problems related with the Heun's equation 
[1]. It is well known in the quantum physics. To the best of our knowledge, this 
phenomenon is observed for the first time in the theory of the perturbations of 
Schwarzschild metric and needs further study to reach the corresponding physical 
understanding. 
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Figure 6: The dependence of Im (0^,2,2) on the mixing angle a G [0, 2n] 



3.3 Perturbations of the Kruscal-Szekeres Manifold 

To study the stable solutions of the first-order Regge- Wheeler perturbation theory 
in this case one has to choose the functions 

<.» =r£>, s > + it>, 8 , 0<i,7(r), 

which are regular at the interior-future-infinity ti nt = 00, and then to impose on 
them the condition T^_(u, s,l) = 0. It is needed to annulate the waves, coming 
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from the space infinity r = oo. Making use of ( 12.31) . (12.51) and (12.71) . it is not hard 
to obtain the following explicit form of this condition: 



lim 

|r|— >oc 



r s+1 e 2i " lar HeunC{2tuj,2s,2iu},2uj 2 , s 2 - 1(1 + l),r) j _ i(5+arau) = 0. (3.6) 



Our numerical study of the case s = I = 2 shows that this equation has no finite 
complex solutions u> with \u>\ < oo. Hence, in this case all existing solutions of RWE 
with no waves coming from special infinity are growing infinitely when approaching 
the point r = 0, i.e. they are unstable in direction of the interior future. 

As an example of such solution one can consider the functions R^~^i( r )- Ac- 
cording to the (12.6bl) . when Im(uj) > these functions are finite and continuous in 
the whole interval r £ (0, oo). They tend continuously to zero at both sides of the 
horizon r = 1 (which in this case is a branching point in C r ) and to infinity - when 
approaching the singular points r = and r = oo. Unfortunately, these solutions 
contain both going to space infinity and coming from there waves. Our numerical 
study of the case s = 1 = 2 shows that it is impossible to fulfill the condition 



(r 



(3.7) 

r = \r\e~ i(K/2+arg(uj)) 



lim r 1+s HeunC(-2iuo, 2iu, 2s, -2a; 2 , s 2 -1(1 + 1) + 2uj 2 , 1 - r) 

\r\—*oo 



^,_|y,| e — i(iv /2 J rarg(uj)) 

= 0, 



which excludes the coming from the space infinity waves and is expecting to yield 
the corresponding spectrum of the frequencies uj. 

Another example are solutions Ruti( r )> considered in the whole interval r G 
(0, oo). Our consideration in Section 3.1 shows that in this case one is able to 
exclude the coming from space infinity waves, obtaining this way the standard QNM 
frequencies. These solutions go to infinity both at the horizon and at the singular 
points r = and r = oo. Hence, they are not stable in direction of the interior 
future and are not suitable for constructing perturbation theory in vicinity of all of 
the three singular points of the RWE. 



3.4 Perturbations of the Exterior of Compact Matter Ob- 
jects 

Here we consider a one-parameter family of real intervals (r*,oo) C (l,oo) with 
Dirichlet's boundary condition at the left end: <& W)8) j(t, r*) = 0. Physically this 
means that the waves are experiencing a total reflection at the spherical surface 
with the area radius r = r* > 1, instead of going freely through it without any 
reflection, as in the case of black holes. 

One can consider this surface as a boundary of some massive spherically sym- 
metric body. Owing to the Birkhoff theorem, in general relativity the gravitational 
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field of any massive body with Keplerian mass M outside the radius r* coincides 
with the corresponding field of a black hole. Hence, the spreading of the waves in 
the outer domain of the massive body is governed by the Regge- Wheeler equation 
( II. II) . too, but the boundary conditions must be different. The corresponding effec- 
tive potential V Sj i(r) in the Regge- Wheeler equation ( II. ip with Dirichlet's boundary 
condition and the corresponding spectrum are shown in FigJT] and FigJHl The spec- 





Figure 7: The potential 
V(r), truncated at r* > 1. 
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Figure 8: The dependence of the 
eigenvalues: n — 0, 1, 2; s — I — 2 
of a compact body on the area-radius 
r # of the totaly reflection of the per- 
turbations, and the first six QNM of 
SBH. 

tral equation, solved numerically to reproduce the FigJHl expresses the requirement 
to have at big distances only going to infinity waves and reads: 



lim | ( 

•I— >+oo 



2i*u ln(|r|e- i ^/ 2 + arg( >'»-l) 



X 



HeunC (- 


■2iu, 


+2iu, 2s, 


-2uo 2 ,2u 2 + s 2 - 


-1(1 + 1), l-\r 


e -i(7r/2+arg(t 


J)) ) 


HeunC (- 


-2iu 


—2iuj, 2s, 


-2cu 2 ,2lu 2 + s 2 - 


-1(1+1), l-\r 


e -i(7r/2+arg(tj))'J 


g2i*o; ln(r«- 


4 )x 


HeunC (- 


-2iuj, +2iu, 2s, 


-2lu 2 ,2lu 2 + s 2 


- 1(1 + 1), 1 


-n 




HeunC (- 


-2iu, —2iu, 2s, 


-2u 2 ,2u 2 + s 2 


-1(1+1), 1- 


-n 



- (3- 



0. 



A more realistic model for real spherically symmetric bodies may be obtained 
considering as a boundary condition a partial reflection and partial penetration of 
the waves through the body's surface. 

Such results may help to clarify the physical nature of the observed very heavy 
and very compact dark objects in the universe. Our consideration gives a unique 
possibility for a direct experimental test of the existence of space-time holes. The 
study of spectra of the corresponding waves, propagating around compact dark 
objects, may give indisputable evidences, whether there is space-time hole inside 
such an invisible object. 
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4 Conclusion 



We have demonstrated that the Heun's functions are an adequate powerful tool for 
description of the Regge- Wheeler linear perturbations of the Schwarzschild gravita- 
tional field not only outside the SBH horizon, but in the interior domain as well. 
They give, too, the possibility for exact treatment of different boundary problems for 
the Regge- Wheeler equation, including ones, which are related with matter spheri- 
cally symmetric compact objects. Using explicit technics we have demonstrated the 
limitations of the Regge- Wheeler perturbation theory of SBH. Since this theory has 
no complete basis of solutions, which are finite everywhere in the range r e (0, oo), 
or at least at finite distances, one can conclude that Regge- Wheeler perturbation 
theory is not adequate for an approximate description of the real properties of the 
full nonlinear problem, which describes small deviations of the Schwarzschild metric. 
Thus we have seen explicitly that this theory may have only a limited range of ap- 
plications. The search for a better type of perturbation theory, or the full nonlinear 
treatment of the small perturbations of the SBH remains an open problem. 
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